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We present the cross-section for the production of a Higgs boson at hadron-colliders at next-to- 
next-to-next-to-leading order (N 3 LO) in perturbative QCD. The calculation is based on a method 
to perform a series expansion of the partonic cross-section around the threshold limit to an arbitrary 
order. We perform this expansion to sufficiently high order to obtain the value of the hadronic cross 
at N 3 LO in the large top-mass limit. For renormalisation and factorisation scales equal to half the 
Higgs mass, the N 3 LO corrections are of the order of +2.2%. The total scale variation at N 3 LO is 
3%, reducing the uncertainty due to missing higher order QCD corrections by a factor of three. 

PACS numbers: 12.38.Bx 


The success of the Large Hadron Collider (LHC) ex¬ 
periments in the exploration and interpretation of phe¬ 
nomena at the TeV scale is due, on the one hand, to 
amazing experimental and technological advances and, 
on the other hand, extraordinary progress in perturbative 
QCD. In particular, the discovery of the Higgs boson [T] 
by the LHC experiments has initiated an era of preci¬ 
sion studies of the properties of the Higgs boson, where 
precise theory predictions for Higgs observables play an 
indispensable role. 

The inclusive gluon-fusion cross-section is a prototypi¬ 
cal example of a theoretical input for the interpretation of 
the experimental observations. It enters not only into the 
extraction of the Higgs-boson couplings from the mea¬ 
surements, but it could also play an important role in 
identifying deviations from the Standard Model predic¬ 
tions in Higgs physics. Unfortunately, the theory predic¬ 
tions for the inclusive cross-section are plagued by signif¬ 
icant theoretical uncertainties. Scale variations at NNLO 
indicate that missing higher order effects are of the order 
of ±9% at LHC energies HE], and the size of this un¬ 
certainty is comparable to the experimental uncertainty 
from the LHC Run 1 [4. 5 . Hence, with a few more years 
of data taking the theoretical uncertainty will be domi¬ 
nant, demanding for an update of the current theoretical 
predictions. 

In this context, a vigorous effort has recently been 
made to compute the inclusive gluon-fusion cross-section 
at next-to-next-to-next-to-leading order (N 3 LO) in per¬ 
turbative QCD. The cross-section at N 3 LO receives con¬ 
tributions from many different building blocks, all of 
which have been computed over the last years, at least 
partially. The three-loop corrections to Higgs production 
in gluon fusion have been obtained in ref. j6], and the cor¬ 
rections from the emission of an additional parton at one 


or two loops were computed in ref. 1M2- In order to ob¬ 
tain a finite result, appropriate ultra-violet and infrared 
counterterms need to be included (12UT4] . While all of 
these contributions had been computed in full generality, 
contributions from the emission of two partons at one 
loop and three partons at tree-level had only been com¬ 
puted in an approximate manner. In particular, for these 
contributions the first two terms in the expansion around 
threshold could be obtained H5HISI, confirming previous 
results for logarithmically enhanced terms in the cross 
section ED and resulting in the complete computation of 
the inclusive gluon-fusion cross-section at N 3 LO in the 
soft-virtual UzHE! and next-to-soft approximations [20] . 
Owing to the universality of soft emissions, the previous 
results have sparked various new results for QCD pro¬ 
cesses at N 3 LO in the soft-virtual approximation [221 . 

Despite this progress, the soft-virtual and next-to-soft 
approximations are insufficient to make reliable predic¬ 
tions for the cross-section, owing to a slow convergence of 
the threshold expansion [20]. In this Letter we close this 
gap, and we present the gluon-fusion Higgs production 
cross-section at N 3 LO in perturbative QCD. We empha¬ 
sise that this is the first ever complete computation of a 
cross-section at N 3 LO at a hadron collider. 

We will describe the main result of our computation in 
this Letter, while a detailed account of the mathematical 
and computational methods will be presented elsewhere. 
Here it suffices to say that we work in the framework of 
reverse-unitarity 1231124] , and we perform a complete re¬ 
duction of the cross-section to master integrals, without 
any approximations. For the double and triple-emission 
contributions at N 3 LO, we can derive differential equa¬ 
tions satisfied by the master integrals [2314251 . which we 
solve as generalised power series around the threshold 
limit. In this way, we obtain at least 37 terms in the 
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threshold expansion of each master integral. An impor¬ 
tant part of our computation has been the evaluation of 
the boundary conditions which are needed for solving the 
differential equations for the master integrals. Many of 
the boundary conditions required in this project had al¬ 
ready been derived in the context of the soft-virtual and 
next-to-soft results [IMS mi- Using similar techniques, 
we have computed the remaining few unknown boundary 
conditions for master integrals which start to be relevant 
only at a high order in the threshold expansion. 

Having at our disposal the complete set of master 
integrals as expansions around the threshold limit, we 
can easily obtain the cross-sections at N 3 LO for all par- 
tonic channels contributing to Higgs production via gluon 
fusion. The partonic cross-sections are related to the 
hadronic cross-section at the LHC through the integral 

* = £/ d,XidX2fi(xi,lJ,f)fj(X2,fjLf)&ij(z,II, r ,fJ,f), ( 1 ) 


where the summation indices i. j run over the parton fla¬ 
vors in the proton, /,; are parton densities and cry are 
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partonic cross-sections. Furthermore, we define z = —j 3 -, 
where m# is the mass of the Higgs boson and y/s is the 
partonic center-of-mass energy, related to the hadronic 
center-of-mass energy VS through s = X\ Xn S. The 
renormalisation and factorisation scales are denoted by 
Hr and iif. We work in an effective theory approach 
where the top-quark is integrated out. The effective La- 
grangian describing the interaction of the Higgs boson 
and the gluons is, 

C eS =~HG%G a »\ (2) 

where H is the Higgs field, G" jJV is the gluon field 
strength tensor and C the Wilson coefficient, known up 
to N 4 LO j2S]. We expand the partonic cross-sections 
into a perturbative series in the strong coupling constant 
evaluated at the scale g r , 


@ij — cr 0 



(3) 


In this expression do denotes the leading order cross- 
section, and the terms through NNLO in the above ex¬ 
pansion have been computed in [231 ffl El]. The main 
result of this Letter is the result for the N 3 LO coeffi¬ 
cient, corresponding to £ = 3 in eq. for all possible 
parton flavours in the initial state. We cast the N 3 LO 
coefficients in the form 
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where we introduce the truncated threshold expansions 
defined by 
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FIG. 1: The N 3 LO correction from the gg channel to the 
hadronic cross-section as a function of the truncation order 
N in the threshold expansion for the scale choice g = mn- 


Here, dgy denotes the soft-virtual cross-section at N 3 LO 
of ref. [171119] and IV = 0 is the next-to-soft approxi¬ 
mation of ref. jtjj- Using our method for the threshold 
expansion of the master integrals, we were able to deter- 
mine the c tJ analytically up to at least n = 30. Note that 
at any given order in the expansion these coefficients are 
polynomials in log(l — z). While this approach does not 
cast the partonic cross-sections in a closed analytic form, 
we argue that it yields the complete result for the value 
of the hadronic cross-section. In Fig. [I] we show the con¬ 
tribution of the partonic cross-section coefficients N 3 LO 
to the hadronic cross-section for a proton-proton collider 
with 13TeV center-of-mass energy as a function of the 
truncation order N. We use NNLO MSTW2008 [21] par- 
ton densities and a value for the strong coupling at the 
mass of the Z-boson of a s (mz ) = 0.117 as initial value 
for the evolution, and we set the factorisation scale to 
g,f = mn- We observe that the threshold expansion sta¬ 
bilises starting from N = 4, leaving a negligible trun¬ 
cation uncertainty for the hadronic cross-section there¬ 
after. We note, though, that we observe a very small, 
but systematic, increase of the expansion in the range 
N G [15,37], as illustrated in Fig. [l] We have observed 
that a similar behaviour is observed for the threshold 
expansion at NNLO. The systematic increase originates 
from values of the partonic cross-section at very small z. 
Indeed, this increase appears only in the contributions 
to the hadronic cross-section integral for values z < 0.1. 
It is natural that the terms of the threshold expansion 
computed here do not furnish a good approximation of 
the hadronic integral in the small z region due to the di¬ 
vergent high energy behaviour of the partonic cross sec¬ 
tions [22]. However, it is observed that this region is 
suppressed in the total hadronic integral and for z < 0.1 
contributes less than 0.4% of the total N 3 LO correction. 
The same region at NLO and NNLO, where analytic ex- 
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FIG. 2: Scale variation of the gluon fusion cross-section at 
all perturbative orders through N 3 LO. 

pressions valid for all regions are known, is similarly sup¬ 
pressed. We therefore believe that the uncertainty of 
our computation for the hadronic cross-section due to 
the truncation of the threshold expansion is negligible 
(less than 0.2%). 

In Fig. [2] we present the hadronic gluon-fusion Higgs 
production cross-section at N 3 LO as a function of a com¬ 
mon renormalisation and factorisation scale /r = /r r = 
///. We observe a significant reduction of the sensitiv¬ 
ity of the cross-section to the scale a 1 . Inside a range 


H € [^.nifl] the cross-section at N 3 LO varies in the 
interval [—2.7%, +0.3%] with respect to the cross-section 
value at the central scale /i = For comparison, we 
note that the corresponding scale variation at NNLO is 
about ±9% [2{ [3j . This improvement in the precision of 
the Higgs cross-section is a major accomplishment due to 
our calculation and will have a strong impact on future 
measurements of Higgs-boson properties. Furthermore, 
even though for the scale choice /i = the N 3 LO cor¬ 
rections change the cross-section by about +2.2%, this 
correction is captured by the scale variation estimate for 
the missing higher order effects of the NNLO result at 
that scale. We illustrate this point in Fig. [3j where we 
present the hadronic cross-section as a function of the 
hadronic center-of-mass energy v / S at the scale /j, = rn ^ L . 
We observe that the N 3 LO scale uncertainty band is in¬ 
cluded within the NNLO band, indicating that the per¬ 
turbative expansion of the hadronic cross-section is con¬ 
vergent. However, we note that for a larger scale choice, 
e.g., p. = mu , the convergence of the perturbative series 
is slower than for At = . 

In table [T] we quote the gluon fusion cross section 
in effective theory at N 3 LO for different LHC energies. 
The perturbative uncertainty is determined by varying 
the common renormalisation and factorisation scale in 
the interval around and in the interval 

around rrin- 


cr/pb 

2 TeV 

7 TeV 

8 TeV 

13 TeV 

14 TeV 

a — vm. 

At = m H 

o qq+ 0.43% 
U,yy -4.65% 

0 04+4.87% 
U - y ^-7.35% 

i c oi +0.31% 
10.0±_3 og % 

14 84 +3 ' 18 2 

19 47+0.32% 

IC/.'H; < —2.99% 

1 & qn+3.08% 

i».yu_ 5 02% 

44 3i+°- 31% 

^±.01 —2 64% 

43 i4+ 2 - 71% 
4 45 % 

49 . 871 ™?$ 

48-57 


TABLE I: The gluon fusion cross-section in picobarn in the effective theory for different collider energies in the interval 
[p^-, rriH ] around n = and in the interval [+r-, 2m h] around /r = . 


Given the substantial reduction of the scale uncertainty 
at N 3 LO, the question naturally arises whether other 
sources of theoretical uncertainty may contribute at a 
similar level. In the remainder of this Letter we briefly 
comment on this issue, leaving a more detailed quantita¬ 
tive study for future work. 

First, we note that given the small size of the N 3 LO 
corrections compared to NNLO, we expect that an esti¬ 
mate for the higher-order corrections at N 4 LO and be¬ 
yond can be obtained from the scale variation uncer¬ 
tainty. Alternatively, partial N 4 LO results can be ob¬ 
tained by means of factorisation theorems for thresh¬ 
old resummation. However, we expect that the insight 
from resummation on the N 4 LO soft contributions is only 


qualitative given the importance of next-to-soft, next-to- 
next-to-soft and purely virtual contributions observed at 
N 3 LO, as seen in Fig. [l] 

Electroweak corrections to Higgs production have been 
calculated through two loops in ref. [32], and estimated 
at three loops in ref. 35j. They furnish a correction of 
less than +5% to the inclusive cross-section. Thus, they 
are not negligible at the level of accuracy indicated by 
the scale variation at N 3 LO and need to be combined 
with our result in the future. Mixed QCD-electroweak 
or purely electroweak corrections of even higher order 
are expected to contribute at the sub-percent level and 
should be negligible. 

Next, we have to comment on our assumption that the 
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FIG. 3: The gluon fusion cross-section at all perturbative or¬ 
ders through N 3 LO in the scale interval [^p,m/r] as a func¬ 
tion of the center-of-mass energy VS. 

top-quark is infinitely heavy and can be integrated out, 
see eq. ([2]). Moreover, we assumed that all other quarks 
have a zero Yukawa coupling. Finite quark mass effects 
are important, but it is sufficient that they are inlcuded 
through NLO or NNLO. Indeed, finite quark-mass effects 
have been computed fully through NLO in QCD 150] , 
while subleading top-quark mass corrections have been 
computed at NNLO systematically as an expansion in 
the inverse top-quark mass [34]. In these references it 
was observed that through NLO finite quark mass ef¬ 
fects amount to about 8% of the K-factor. At NNLO, 
the known —— corrections affect the cross-section at 

^top 

the ~ 1% level. A potentially significant contribution 
at NNLO which has not yet been computed in the lit¬ 
erature originates from diagrams with both a top and 
bottom quark Yukawa coupling. Assuming a similar per¬ 
turbative pattern as for top-quark only diagrams in the 
effective theory, eq. ([2]), higher-order effects could be of 
the order of 2%. We thus conclude that the computation 
of the top-bottom interference through NNLO is highly 
desired in the near future. 

Finally, the computation of the hadronic cross-section 
relies crucially on the knowledge of the strong coupling 
constant and the parton densities. After our calculation, 
the uncertainty coming from these quantities has become 
dominant. Further progress in the determination of par- 
ton densities must be anticipated in the next few years 
due to the inclusion of LHC data in the global fits and the 
impressive advances in NNLO computations, improving 
the theoretical accuracy of many standard candle pro¬ 
cesses. 


To conclude, we have presented in this Letter the 
computation of the gluon-fusion Higgs production cross- 
section through N 3 LO in perturbative QCD. While a 
thorough study of the impact of electroweak and quark 
mass effects is left for future work, we expect that the re¬ 
maining theoretical uncertainty on the inclusive Higgs 
production cross-section is expected to be reduced to 
roughly half, which will bring important benefits in the 
study of the properties of the Higgs boson at the LHC 
Run 2. Besides its direct phenomenological impact, we 
believe that our result is also a major advance in our un¬ 
derstanding of perturbative QCD, as it opens the door to 
push the theoretical predictions for large classes of inclu¬ 
sive processes to N 3 LO accuracy, like Drell-Yan produc¬ 
tion, associated Higgs production and Higgs production 
via bottom fusion. Moreover, on the more technical side, 
our result constitutes the first independent validation of 
the gluon splitting function at NNLO [14], because the 
latter is required to cancel all the infrared poles in the 
inclusive cross-section. In addition, we expect that the 
techniques developed throughout this work are not re¬ 
stricted to inclusive cross-sections, but it should be pos¬ 
sible to extend them to certain classes of differential dis¬ 
tributions, like rapidity distributions for Drell-Yan and 
Higgs production, thereby paving the way to a new era 
of precision QCD. 
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